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ABSTRACT 

Near-UV observations of the planet host star WASP- 12 uncovered the apparent 
absence of the normally conspicuous core emission of the Mg Ilh&k resonance lines. 
This anomaly could be due either to (1) a lack of stellar activity, which would be 
unprecedented for a solar- like star of the imputed age of WASP- 12; or (2) extrinsic ab- 
sorption, from the intervening interstellar medium (ISM) or from material within the 
WASP- 12 system itself, presumably ablated from the extreme hot Jupiter WASP- 12 b. 
HIRES archival spectra of the Ca II H&K lines of WASP- 1 2 show broad depressions in 
the line cores, deeper than those of other inactive and similarly distant stars and similar 
to WASP-12's Mg Ilh&k line profiles. We took high resolution ESPaDOnS and FIES 
spectra of three early-type stars within 20' of WASP- 12 and at similar distances, which 
show the ISM column is insufficient to produce the broad Ca II depression observed 
in WASP- 12. The EBHIS H I column density map supports and strengthens this con- 
clusion. Extrinsic absorption by material local to the WASP- 12 system is therefore the 
most likely cause of the line core anomalies. Gas escaping from the heavily irradiated 
planet could form a stable and thick circumstellar disk/cloud. The anomalously low 
stellar activity index (\ogR HK ) of WASP- 12 is evidently a direct consequence of the 
extra core absorption, so similar HK index deficiencies might signal the presence of 
translucent circumstellar gas around other stars hosting evaporating planets. 

Subject headings: stars: individual (WASP-12) — ISM: clouds — Planets and satellites: 
individual (WASP-12) — Planet-star interactions 
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1. Introduction 



WASP- 12, a late F-type main sequence star younger than 2.65 Gyr (|Fossati et al.ll2010b|) . 
hosts one of the hotte st, most bloated known exoplanets, orbiting extremely close to the host 



star (HebbetaL 



2009). WASP-12b is one of the most irradiated transiting planets known, and 



provides key tests of the "blow-off" hypothesis whereby a significant fraction of the planet's 
atmosphere is being eroded and escapes. 



Fossati et al. 



(2010a) and 



Haswell et al 



(120121) analysed Hubble Space Telescope (HST) 



near-UV (NUV) spectra obtained with the Cosmi c Origin Spectrograph (COS), confirming th e 



on-going blow-off, similar to that of HD209458 (IVidal-Madjar et al. 



2003; 



Linskv et al. 



2010) 



NUV dips during the planetary crossings are up to three times deeper than the optical transits 
of WASP- 12 b , suggesting there is diffuse gas extending well beyond the planetary Roche lobe. 



Li et al 



(2010) suggested that the gas lost by the planet could form an accretion disk around the 
star. The HST data revealed a variable and early ingress, contrary to model expectations, but with 
egresses compatible wit h the optical ephemeris. The early- i ngress phenomenon has been modeled 



by various authors (e.g., 



Lai et al. 



20101: 



Vidotto et al. 



2010 



Bisikalo et al. 



do not yet explain the presence of an ingress as early as orbital phase -0.17 (|Haswell et al. 



2013), but these mode 



2012). 



The HST data further disclosed a remarkable anomaly in the stellar NUV spectrum: the 
apparent absence of the normally bright emission cores in the Mg II h&k resonance lines 



(Haswell et al. 



2012|) . Significantly, the anomaly is always present, regardless of the planet's 
orbital phase. The absence of Mg II c ore emission was completely unexpected g iven the spectra l 



type and imputed age of the star (e.g., 



Guinan & Engle 



2009). Figure 16 of 



Haswell et al. 



(2012) 



compares the near-UV spectra of WASP-12, HD102634, HD107213 and Procyon; Procyon is a 
mid-F subgiant commonly used as a n inactive comparison s tar, while HD102634 and HD107213 
are similar in temperature and age (|Valenti & Fisherll2005|) to WASP-12. The comparison 
highlights the complete lack of emission in the WASP-12 line cores. 



Fossati et al.l ((2010b) and IHaswell et al.l (|2012f) discussed scenanos that could lead to the 
anomalous line cores: (1) intrinsically low stellar activity; or (2) extrinsic absorption, either from 
the i nterstellar med i um ( ISM) or from material within the system, for example in a circumstellar 



disk. 



Has well et al 



(|2012D showed that the ISM column density required to reproduce the apparent 
Mg II absorption was 10 times greater than expected along the WASP- 12 line of sight, but that 
planetary mass loss potentially could produce the large column required to extinct the Mg II cores. 

Opportunities to use the Mg Ilh&k lines to compare WASP- 12 with other stars are few: 
inactive distant stars have rarely been observed at high resolution in the NUV, because of the 
bias of interest in NUV-bright stars, which either are active or nearby. The optical Ca IIH&K 
resonance lines also display core em ission of ch romospheric origin closely related to that of the 



2008). The Ca II lines have been used to assess the 



Mg II lines, although not as strong (Hall 
activity of extensive samples of late-type stars, including planet-hosts, some of which indeed are 



inactive and distant (Knutson et al 



2010). 



In Sect.[2]we compare WASP-12's Ca IIH&K line profiles with those of other inactive and 
distant planet-hosting stars; in Sect. [3] we evaluate the ISM absorption along the WASP- 12 line 
of sight; and in Sect. @]we discuss possible origins of the anomalously depressed resonance line 
cores. 



2. The Ca II H&K line profiles: activity and anomaly 

It is well known that chromospheric activity correlates with the stellar rotational velocity. 
For WASP- 12 the rotation rate is unknown, but the projected rotational velocity (v sin i) has an 
upper limit of only 1 km s _1 (E. Simpson, private communication). Because a low v sin i could 
potentially imply an intrinsically low level of stellar activity, we must consider whether depressed 
chromospheric emission might be responsible for the Mg II core anomaly. Ideally, we would 
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compare with stars of similar type which do not host hot-Jupiters. However, planet-hosts are 
far better studied at high spectral resolution in the visible, so we exploited planet-host spectra 
obtained for planet detection through Doppler-reflex or stellar characterisation. 

Figured] compares the WASP-12 Ca IIH&K lines with those of HAT-P-7 and WASP-1. The 
spectra were acquired with the High Resolution Echelle Spectrometer (HIRES) on the KECK I 
telescope. We obtained the reduced dataf from the HIRES archive;;, and focused on spectra 
with the highest signal-to-noise ratio (S/N) and taken with the same instrument settings as the 
out-of-transit phases of WASP-12. 

The fundamental parameters of the three stars are listed in Table [T] Within the uncertainties, 
the comparison stars are comparable in effective temperature (Teg) to WASP-12. Reliable ages 



are not available for HAT-P-7 and WASP- 1 , but their low v sin 



j and \ogl( HK values indicate they 



certainly are not young. Distances were estimated based on the V-ba nd magnitudes, r p ff ; and r adii 



reddenings E(B - V) then were derived from the extinction maps of lAmdres & Lepind (|2005[) 



'htt p : //www . astro . caltech . edu/~tb/makee/| 
2 http: //www2 .keck.hawaii . edu/koa/public/koa.php 

3 The v sin i values reported in Table \T\ should be considered upper limits as 
ignored possible contributions from "macroturbulence". 



Torres et al 



(2012) 
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Fig. 1. — Comparison between the Ca IIH&K line profiles (right and left panels, respectively) 
of WASP-12 and HAT-P-7 (upper panels) or WASP-12 and WASP-1 (lower panels). Horizontal 
bars indicate spectral windows used to normalise the flux profiles. Short vertical lines in the right 
panels mark the location of hydrogen He. The per pixel uncertainties are of the order of 3%. The 
lower panel within each plot shows the ratio of the spectra with WASP-12's as the numerator. In 
each case WASP-12 has extra absorption throughout the inner ~4 A of the line profile. 
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Because WASP- 12 has the lowest vs'mi value, we broadened its spectrum to match the 
comparison tracings. Flux normalisation was performed in two steps. We removed the blaze 
function adjusting the continuum using a synthetic spectrum with atmospheric parameters 
appropriate to the three stars. We then matched the fluxes of the far line wings to WASP- 12, for 
the H and K lines separately, using the windows indicated in Fig.Q] The latter adjustments were 
all <5%. 



As anticipated from the earlier Mg Ilh&k line observations (|Haswell et al.ll2012|) . the 
Ca II H&K line cores of WASP- 1 2 are significantly deeper than those of HAT-P-7 and WASP- 1 , 
despite similarities in the stellar properties. The anomaly covers a ~4 A wide region, similar 
to the Mg Ilh&k depression. The same comparison relative to CoRoT-1, HAT-P-8, HAT-P-16, 
HAT-P-24, HAT-P-33, HAT-P-39, Kepler-8, Kepler-25, TrES-4, and WASP-3, always shows the 
central ~4 A of the line core the most depressed in WASP- 12. The flux adjustment might seem 
to be challenging in this heavily blanketed spectral region, but there are enough relatively clean 
intervals in the Ca II wings that our results are very little affected by the normalization procedure. 
We carefully checked that background light under-/over-correction could not explain the observed 
anomaly in the WASP- 12 line core. 



Knutson et al 



(|2010f) measured the logi?^ activity index in 50 transiting planet host 



stars, finding that WASP- 12 was the least active, by that measure, with a remarkably low 
logR' HK = -5.500. Figure |2] shows that WASP- 12 is an extreme out: 



activity (log R HK ) 



2008, 



20111) . and 



plane populated by t he samples from 



Knutson et al. 



Wright et al. 



(2004), 



color (B - V) 


vs. 


Jenkins et al. 


(2006 



(120101) . As the logR HK index is logarithmic, WASP-12 indeed 



falls well below the lower boundary of the envelope of activit y displayed by the stars con sidered 



here. For solar metallicity dwarfs the lower bound is at -5.1 (|Henry et al 



1996; 



Wright 



2004) 



and it decre ases slightly to -5.15 for stars with the somewhat enhanced metallicity of WASP-12 



(Saar 



201 lb . The lower bound of the index results from a "basal" level of chromospheric emission 
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which seems to depend mainly on T e s and only weakly on metallicity (as noted above). Stars 
with indices above the basal level are affected by varying degrees of dynamo magnetic activity, 
which does correlate with the rotational velocity. The basal flux level should be a hard lower 
bound on the chromospheric emission of a late-type main sequence star, and any otherwise normal 
dwarf that falls below the basal level must be affected by extrinsic absorption which attenuates 
the intrinsic core emission. To be sure, evolved subgiants can fall below the basal boundary, as 



discussed by 



Wrightl (|2004|) . Figure |2] contains contamination by such subgiants, but they can 



be di scriminated by the broadening of their Ca II emission cores by the "Wilson-Bappu effect" 



(e.g. 



Avres 



1979|). In short, as suggested in 



low \ogR HK index of WASP- 12 derived by 



Haswell et al 



Knutson et al 



(2012), it is likely that the abnormally 



(|2010l) is simply biased by extrinsic 



absorption, hiding the chromospheric signatures that must be present intrinsically in Ca II H&K. 
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Table 1: Fundamental parameters of WASP- 12 and two comparison stars. 



Star 


T e ff 


v sin i 


[Fe/H] 


logR HK 


distance 


E(B - V) 


name 


[K] 


[kms" 1 ] 


[dex] 




[pc] 


[mag] 


HAT-P-7 


6441 


5.0 


+0.15 


-5.018 


310 


0.027 


WASP-1 


6160 


1.7 


+0.14 


-5.114 


388 


0.020 


WASP- 12 


6250 


<1 


+0.28 


-5.500 


439 


0.144 



Not e. — r e ff, vsini, and metallicity values are from lTorres et al.l (120121) ; logR HK values are from iKnutson et al. 
(2010); distances and reddening were derived as explained in the text. 
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3. ISM measurements along the WASP- 12 sightline 



To determine the contribution of interstellar absorption to the depressions in the cores of 
both Mg II and Ca II, we carried out a s tudy of the ISM column in the WASP- 12 sightline. 



We used ESPaDOnS (Donati et al. 



19971) at the Canada-France-Hawaii Telescope (CFHT) and 
FIES at the Nordic Optical Telescope (NOT) to obtain high resolution spectra of three early-type 
stars (HD257926, HD258049, HD258439) lying within 20' of WASP-12, and at about the same 
distance. We observed hot stars for three reasons: 1) they are brighter than WASP-12 at the 
same distance, allowing high-S/N measurements; 2) they are typically fast rotators, simplifying 
identification of the sharp ISM absorptions; and 3) they do not exhibit stellar activity, which might 
otherwise affect interpretations of the Ca II H&K lines. 

The ESPaDOnS spectra cover 3700-10400 A, with a resolving power of 80000 in the "star 
only" instrumental configuration, while the FIES spectra cover 3700-7400 A, with a resolving 
power of 67 000. The spectra revealed that HD258049 is a spectroscopic binary, while HD258439 
is a chemically peculiar star (possibly magnetic), and therefore less straightforward to model. For 
this reason we focused on HD257926, apparently a chemically normal early A-type star. It lies 
15.7' away from WASP-12 on the sky, at about twice the distance (~650pc). We estimated the 
fundamental parameters by comparing observed hydrogen line profiles (Ha, H/3, Hy, and US) with 
synthetic spe ctra calculated with the SYNTH3 code (|Kochukhovll2007|) . utilizing LLmodels model 



atmospheres (IShulyak et al 



20041) . We obtained: r eff ~10700K, logg~4.25, vsini -120 km s" 1 . 



Narrow ISM absorption lines are clearly visible in HD257926's Ca II H&K and Na ID1 
and D2 lines: see Fig. [3] (top-left panel) for Ca II K. Each doublet member shows at least two 
ISM velocity components. For conciseness we focus on the higher S/N ESPaDOnS spectrum, 
but the FIES spectrum provided comparable results. We removed the stellar contribution using a 
synthetic spectrum, and fitted multiple Gaussians to the residual ISM absorptions. Table[2]lists the 
equivalent widths of the ISM lines. We calculated synthetic ISM line profiles for a range of "b" 
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Fig. 2.— WASP- 12 (blue triangle ) in B - V vs. l ogi?' plane. 



V < 1 .0 observed by 



Kn 



by 



utson et al 



Wright et al 



(12004 (squares), 



Jenkins et al. 



compared with sta rs with B - 



(2006 



2008 



20111) (stars), and 



(2010) (red pluses). Color and activity indices adopted for WASP- 12 are those given 



Knutson et al.1 (|2010h . Circles indicate positions of planet hosting stars XO-4, CoRoT-1, WASP- 



13, WAS P- 1 7 , and WASP- 1 8 (see Sect. U). The dotted line indicates the minimu m activity valu e 



within the 



Wright et al. 



(2004) sample, accounting for contamination by subgiants (Wright 



2004). 
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parameters (characterizing thermal and turbulent line broadening) and column densities, and for 
each b value determined the column density that reproduced the observed equivalent width. For 
Ca II, we then compared the derived H and K column densities, and determined the b value for 
which the two column densities agreed. We proceeded in the same way for the Na I D lines. In all 
cases we allowed the two Gaussian profiles to have variable widths. Our conclusions were not 
significantly affected by these choices. 
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Uncertainties were based on (S/N) 2 values of the ESPaDOnS spectrum of HD257926, 
separately in the regions of the Ca II and Na I lines. We then applied a Monte Carlo scheme to 
estimate the parameter uncertainties in the Gaussian fitting. 

In estimating the final Ca II and Na I column densities, we considered the effect of the 
equivalent width uncertainties and the influence of the b parameter, which was varied in the 
2. 8-4.0 km s" 1 rang e, compatible with the spread of values observed in the direction of WASP- 12 



(IWelsh et al. 



2010J). For the Ca II lines, we derived b values of ~3.5 km s _1 , while for the Na I 
lines we derived b values of ~4 and <3 km s" 1 for the weak and strong components, respectively. 
The smaller b value of the stronger Na I components suggests a lower temperature in the ISM 
clouds producing these absorptions. 

Table [2] summarises the final results. The colum n densities and b parameters derived for 



Welsh et al 



(2010) for stars at a distance of 



Ca II and Na I are consistent with those reported by 
~400pc in a direction close to that of both HD257926 and WASP- 12. The firm conclusion, as 
illustrated graphically in the lower panel of Fig. |3] is that the ISM absorption along the WASP- 12 
sightline is insufficient to cause the apparent broad absorption seen in the Mg II and Ca II cores. 
Figure |3] also shows the position of a narrow absorption feature in WASP-12's line profile, which 
is absent from the ISM absorption. It is blue-shifted by 14.5 and 14.2 km s -1 for the Ca IIH&K 
lines, respectively. This feature is absent in all the other stars we examined in Sect. EI 

Nevertheless ISM absorption might vary over small angu lar distance ? but the total H I colum n 
density maps from the Effelsberg-Bonn HI Survey (EBHIS; 



Winkel et al. 


2010: 


Kerp et al. 


2011) 



show (Fig. O that the H I column density variation within 20' from WASP- 12 is small and 
WASP- 12 has a lower column than HD257926. 
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Fig. 3. — Top-left panel: HD257926 (solid line), observed with ESPaDOnS; synthetic spectrum 
(dash-dotted line) used to recover the ISM Ca II K absorption profile (dashed line). Bottom-left 
panel: WASP- 12 (solid line) compared to the Sun (dash-dotted line: National Solar Observatory 
Digital Archive). The Solar spectrum was normalised with the procedure of Fig. [TJ and convolved 
to match the HIRES resolution. Note the excellent agreement in the far wings. The dashed line 
applies ISM absorption obtained from HD257926 to the Solar spectrum. The dashed vertical line 
indicates the position of the ISM-unrelated feature in WASP-12's spectrum. Right panel: EBHIS 
total H I column density map in the WASP-12's direction. Cross and plus mark the positions of 
WASP- 12 and HD257926, respectively. 
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Table [2] also lists column densities obtain ed assuming average interstella r conditions, 



N ff =5 .8xl0 21 xElg-V") dSavage & Mathis 



1979), the 



(2003), and either the gas-phase ISM abundances of 



abundances of 



SM ionization mixi n g by 



Asplund et al. 



Savage & Sembachl dl996|) or the solar 



'risch & Slavin 



(120091) . The measured Na and Ca ISM column densities in the 
direction of WASP- 12 apparently are slightly larger than t he average ISM, comp atible with the 



presence of a relatively dense H I region in that direction (IBen Bekhti et al. 



2012). We estimated 



the corresponding Mg II column density assuming the solar a bundance, obtaining log N= 16.42 



in agreement with the value proposed by 



Haswell et al. 



(201 21), but far below the logN=17.30 



required to reproduce the full Mg II line core depressions (Ha swell et al 



2012 ) 



4. Discussion 



In Sections [2] and [3] we showed that neither lack of activity nor strong ISM a bsorption is 



ikely to be responsible for the anomalously deep resonance line cores in WASP- 12. 



Haswell et al 



(120 12|) suggested that the anomaly mig ht be caused by extri nsic absorption by diffuse material 



within the WASP- 12 planetary system. 



Fossati et al. 



(12010b|) searched for the infrared signature 



of a cool component to WASP-12's SED, but found none. 



It is we 



Fossati et al 



1 estab l ished that WASP- 12 b is experiencing a blow-off phase (e.g. LLi et al]|2010l : 



2010al : 



Haswell et al 



2012J). Thus, the planetary atmosphere would be an obvious 



reservoir of gas to fill a circumstellar disk/cloud. Given the variability of chromospheric activity 
and that the planet's orbit is only one stellar diameter away from the surface of WASP- 12, there 
are many processes stripping gas from the planet, and subsequently acting upon it, making the 



geome try complex and time variable. The broadband NUV transits reported by 



Haswell et al. 



(|2012l) directly revealed dramatic variability in the diffuse gas distribution. 



Li et al. 



(120 10|) suggested that gas lost by the planet would fall toward the star, forming an 
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Table 2: ISM properties derived from HD257926 





Ca II K 


Ca II H 


NalDl 


Na ID2 


Vi—blue 


+ 1.8+0.1 


+2.0+0.2 


+4.2+0.2 


+4.2+0.2 


Ur—red 


+ 17.3+0.2 


+ 18.0+0.2 


+ 17.0+0.1 


+16.9+0.1 


FWH1VW 


10.0+0.2 


9.8+0.5 


8.0+0.5 


10.5+0.5 


FWHM red 


11.7+0.3 


10.4+0.7 


9.3+0.2 


10.1+0.2 


EQWblue 


82+2 


56+2 


35+2 


79+4 


EQW red 


79+2 


48+3 


182+2 


212+3 


logNite 


12.2+0.1 


11.6+0.1 


logN red 


12.1+0.1 


12.8+0.1 


logN ss96 


9.79 


11.37 


logN s „„ 


13.52 


12.25 



Note. — Radial velocity (f, - in km s _1 ), full width at half maximum (FWHM - in km s~'), equivalent width (EQW 
- in mA), and column density (in cm 2 ) of the Ca II H&K and Na I D 1 and D2 ISM velocity components ("blue" and 
"red"), measured in HD257926. The final two rows list total Ca II and Na I column densities derived for average ISM 
conditions and the WASP-12 redd ening, assuming eithe r the gas-phase ISM abundances of Savage & Sembach (1996 



SS96) or the solar abundances of 



Asplund et al 



(2009) (see text). 
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accretion disk interior to the planet. In contrast, IVidotto et al.1 (|2010h interpret the observations 
in terms of bow-shocked material entrained in the planet's magneto sphere. It is certain that the 
stellar wind and the stellar radiation pressure, along with the dynamic magnetic fields of both 
planet and star, play important roles in determining the time-dependent behaviour of the gas. The 
diffuse gas, which shares the planet's specific angular momentum, but has a greatly enhanced 
spec ific cross-section to radiation, is likely to move outwards from the planet's orbit as suggested 



in 



Haswell et al. 



(|2012|) . In that scenario, it is easy to imagine that the gas will not be tightly 



constrained to the orbital plane, but will disperse to cover mu ch of the dis k of the star from our 



near-edge-on viewing angle. In the accretion disk scenario of 



Li et al 



(2010) it is less obvious that 



the diffuse gas will s hadow the wide spread bright chromospheric "network" on the stellar disk, 



but the simulation by 



Bisikalo et al. 



(|2013|) . suggests that the planet could expel gas to distances 
as large as five planet radii perpendicular to the planet's orbit (about half the stellar diameter, if 
the ejection is symmetrical about the orbital plane). The relatively high temperature of the disk 
should prevent it settling onto the orbital plane. But, if this indeed were to occur, the large mass of 
the material comprising the disk would force the settling time scale to be longer than the orbital 
period, therefore the planet could continually pump up the disk thickness by feeding new material 
into it. 

We have argued that the flux depression anomaly we detected in the cores of the Mg II and 
Ca II resonance lines of WASP- 12 is evidence for the presence of a translucent circumstellar 
disk/cloud of material lost by the evaporating exoplanet. Figure [2] (bottom-left corner and 
insert) highlights five other stars, also hosting hot-Iupiters, that have anomalously low stellar 
activity indices: XO-4, CoRoT-1, WASP- 13, WASP- 17, and WASP- 18. The latter is particularly 
noteworthy given the many similarities with the WASP- 12 system. The only significant differences 
are the planet mass (WASP- 18 b is ~ 7 times more mass ive than WASP- 12 b) and the stellar 



age (WASP- 18 is only -0.6 Gyr old; 



Brown et al. 



201 lh . The heavier companion and larger 



stellar activity would undoubtedly have an impact on the geometry and stability of a planet-fed 
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disk/cloud. A practical advantage is that WASP-18 is 2.3 magnitudes brighter than WASP-12, 
allowing more precise measurements of the planetary exospheric radius, as well as the prospect of 
acquiring far-UV spectra (which would conclusively settle the "activity" issue). Our findings and 
these considerations make the WASP-18 planetary system a key target for future UV observations 
to continue to build the legacy of HST studies of evaporating hot-Jupiters. 
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